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I. INTRODUCTION 


The knowledge of boiling crisis in the nuclear reactors 
during exponential heat supply is important for the safety 
and efficient performance. The moment of critical heat 
flux (CHF) after which sharp reduction in heat transfer 
coefficient (HTC) is observed may lead to the rapid surge 
in the core temperature. Thus the formation of vapor 
blanket at CHF may lead to core meltdown accident. 
Hence understanding of the mechanism of transient CHF 
during exponential power escalation is of paramount 
importance. Researchers have contributed to explaining 
the CHF mechanisms by various approach viz. (i) Kelvin— 
Helmholtz instability between the upward flow vapor 
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Abstract— In present study, the experiments are carried on sample of 
different surface roughness to investigate the transient heat transfer 
Phenomenon at the saturated condition of the distilled water. The surface 
roughness (Ra) ranges from 0.106 um to 4.03 um. The boiling crisis is 
observed during each transient heat supply. The high-speed camera of 
1000 fps is used to observe the stages of boiling during different transient 
and to confirm the moment of critical heat flux (CHF). The empirical 
relation is presented for transient CHF and corresponding heat transfer 
coefficient (HTC). It is found that transient CHF is a function of both Ra 
and y. The hydrodynamic model is developed for prediction of CHF at 
different rate of exponential heat supply for the wide range of Ra by 
incorporating y. 


columns and downward flow liquid, (ii) dry-out of the 
liquid layer i.e. micro/macrolayer dry-out. Kutateladze [1] 
and Zuber [2] considered the Kelvin-Helmholtz instability 
as surface-fluid interaction fails due to relative motion 
between vapor column and surrounding liquid. Chang [3] 
considered the forces acting on the bubble during lift off 
from the horizontal surface and claimed that the critical 
velocity of the upward moving bubble is responsible for 
CHF. Haramura and Katto [4] stated near field evaporation 
phenomena through the macrolayer evaporation model and 
estimated CHF at the dry-out condition of macrolayer. 


Lay and Dhir [5] and Pasamehmetoglu et al. [6] 
developed the theoretical model based on microlayer 
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evaporation. Zhao et al. [7, 8] developed steady-state and 
transient CHF model by considering various boiling 
aspects like time variant microlayer thickness, microlayer 
evaporation, macrolayer evaporation and transient 
conduction over the entire period of boiling. Hence, 
hydrodynamic models were coupled with surface-fluid 
properties like contact angle [9], surface roughness [10,11] 
and heater orientation [12]. The effects of boiling surface 
properties must be included in a robust and widely 
applicable CHF model. Kandlikar [11] proposed a force 
balance model which includes the contact angle of the 
bubble. Ahn et al. [13] included the term for capillary 
wicking in the Kandlikar’s model to estimate the CHF. 
Quan et al. [14] proposed the force balance model, which 
includes roughness factor (r), to predict CHF for the micro 
structured surface. Kim et al. [15] developed the CHF 
model by considering capillary force, surface roughness 
(Ra) and static contact angle. The capillary force through 
the unidirectional scratches was predicted by assuming 
number of capillary tubes underneath the growing bubble. 
It is noticed from the literature that theoretical model can 
be developed from force balance approach to predict the 
CHF. Far-field and near-field models are developed by 
considering surface roughness (Ra), roughness parameter r, 
static contact angle, surface wetting property, surface 
orientation and heater size. In the present study, influence 
of surface roughness R, and period of exponential heat 
supply on transient CHF is studied. A pioneering study is 
carried to include new term, so-called time constant y in 
the CHF model. 


Il. EXPERIMENTATION 
2.1 Experimental Setup 


A boiling chamber with the test section and condenser 
assembly include in the experimental setup. The schematic 
of experimental setup and visualization unit is shown in 
Fig. 1. The detachable top and bottom flange are provided 
to the boiling chamber. The condenser coil is attached to 
the top flange whereas bottom flange can accommodate 
the test section assembly. The bulk fluid temperature (T)) 
and chamber pressure are measured by the thermocouple 
and pressure transducer, respectively. The transparent 
toughened borosilicate glass watch windows of 115 mm 
diameter and 15 mm thickness are provided to the wall of 
boiling chamber to conduct the visualization study by high 
speed camera. The saturation condition of the distilled 
water is maintained by the two high density cartridge 
heaters each of 1000W capacity. The setup is synchronized 
with high speed camera, NI-9213 temperature module and 
NI-9264. 
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Fig. 1: Experimental Setup 


2.2 Test Section 


Thick copper sample of 20 mm length and the diameter of 
20 mm is prepared as shown in Fig. 2. The 840 W high 
density cartridge heater is used to supply the heat. The 
thermal paste is applied on both the contact surfaces and 
thereafter, sample is screwed with the heating block which 
ensures the perfect surface contact between them. The 
glass wool insulation is provided over the sample and 
heating block. An O-ring and high-temperature non- 
corrosive RTV silicone gasket ensures the leak proof 
assembly. Three K-type sheathed thermocouples of 1 mm 
diameter are implanted from the top of the sample surface 
at 2 mm, 6 mm, and 10 mm. 

__ Boiling surface 
4 4 2 j Set —— Gasket 

- Taa 


m 


Tasi 
__ Test Sample 


Threaded Joint 


Cartridge Heater 


__ Heating Block 
Glass wool 


Teflon 


Fig. 2 : Test Section 


2.3 Experimental Procedure 


The assembly of test section is considered as the 
axisymmetric system. The uniform heat flux from the 
boiling surface is assumed due to uniform surface 
characteristics. The heat flux dissipated to the boiling fluid 
and the surface temperature can be measured using 
thermocouple. 


The time variant heat flux from the boiling surface due to 
exponentially varying heat supply is calculated by 
Equation 1. 
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ie 2Ax (1) 
where Ax is the gap between two adjacent implanted 
thermocouples. 


The temperature of the boiling surface is calculated by 
using Equation 2. 

t _ pt —øe (ma 
Te Ta- as ( Kew ) (2) 


where, Xm-1 is the gap between boiling surface and the 
thermocouple (Tin-1), as given in Fig. 2. 


Heat transfer coefficient (HTC) between the boiling 
surface and water is estimated by Equation 3. 


h,, = 3 
ts T-T) ( ) 


The values of uncertainty in the measured and estimated 
parameters are given in Table 1. 


The values of uncertainty in the measured and estimated 
parameters are given in Table 1. 


Table I : Uncertainties of measured and calculated 


parameters 
Parameter Uncertainty 
T +0.1°C 
x +0.0001 m 
P 0.2 % 
heat flux 11.91 % 
HTC 16.20 % 


Il. RESULTS AND DISCUSSION 


The roughness parameters of the test sample, which were 
measured before and after boiling tests, are tabulated in 
Table 2. The Ra value of the test sample is considered as 
the roughness parameter in this study. 


Table 2 : Roughness parameters in um 


Ra R: R Sm 
0.106 1.20 0.14 13.2 
0.83 7.05 1.07 26.8 
1.87 13.30 2.40 35.7 
3.17 22.91 4.12 42.2 
3.59 27.55 4.09 44.8 
4.03 26.50 4.95 45.2 


The experiments are conducted on different samples with 
exponentially varying heat supply and respective boiling 
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curves are obtained. Effect of roughness on transient 
boiling heat transfer is noticed in the boiling curves 
obtained for y=3, as shown in Fig. 3. It is clearly found that 
boiling curves moves onto the left with increase in values 
of Ra, indicating that the heat transfer enhancement is due 
to the rough surface. The heat transfer enhancement due to 
rough surface can be justified by the mechanism of liquid 
replenishment and surface wettability. Unidirectional 
scratches made on the surface act as a passage for liquid 
supply to the nucleation sites. The wetting property of the 
surface also improved due to capillary wicking along the 
scratches. The irregularities increase with Ra which has 
more potential to be nucleation sites. Thus increased 
number of bubble and the improved liquid supply 
mechanism to the nucleation sites resulted in heat transfer 
augmentation. Also, the nucleation boiling temperature is 
found to be decreased with increase in values of Ra. The 
activation of pre-existing vapor in the cavities of rough 
surface results in formation of bubble at lower wall 
superheat. 


— R,=0.106 ym 
=- - R,=0.83 um 
++ R =1.87 pm 
=) = R,=3.17 um 
== R,=3.59 pm 


q, (kW/m°) 


Fig. 3: Boiling curves of the test samples at y=3 


It is observed during visualization study that exponentially 
increasing heat supply creates the instability in the bubble 
formation. The waiting period between two successive 
bubbles drops drastically with increase in rate of heating. 
Thus, sudden increase in the bubble frequency results in 
the vertical bubble coalescence. Formation of vapor 
column is observed during visualization at this stage, as 
shown in Fig. 4. This mechanism of bubble formation 
decreases the heat transfer rate from the surface. 
Hydrodynamic instability affects the liquid replenishment 
adversely which turn to horizontal bubble coalescence. 
Thus, the blanket of vapor occupies the entire boiling 
surface. The transition mechanism of nucleated bubble 
from the fully developed nucleate boiling (FDNB) regime 
to the film boiling results in drop of HTC. This transition 
point of boiling curve which corresponds to the maximum 
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HTC is identified as CHF. The transient CHF and 
corresponding HTC is presented in Fig. 5 and Fig. 6, 
respectively, at different Ra and y. The values of g curus 
are plotted against Ra and y separately and exponents are 
obtained. The final relation for q’cur.s as a function of Ra 
and y is derived as given in Equation 4 by least square 
regression. It is believed that maxs is a function of q` cris, 
Ra and y. The relation developed for Amaxis is given in 


Equation 5. 
qdZar se = 1031.9 (y) °**(R,)*™* o 
S 104 
__ Sane = 0.04 (y) aia al OAT (5) 


t=3min 04sec) | t=3min 25sec ' t=3min 42sec 


Onset of boiling Nucleate boiling regime 
t=4min 05sec t=4min 35sec t=Smin 10sec 


A 


Vertical coalescence 
t=8min 42sec! 


Fig. 5 Variation in q cur«s with Ra and y 
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Fig. 6 Variation in hmaxts with Ra and y 


IV. HYDRODYNAMIC MODEL 


The bubble formed on the upward facing rough surface 
experiences the forces as shown in Fig. 7. Force due to 
change in the momentum (FM) acts on the bubble 
circumference which pulls the bubble along the surface. 
Simultaneously, bubble position retains due to forces like 
surface tension force, gravity and capillary force acting on 
it. It is observed during present study that horizontal 
coalescence is responsible to turn FDNB to film boiling 
regime. This transition in phase of boiling occurs when 
FM surpasses the sum of drag forces. The moment of CHF 
is identified by the force balance as given in Equation (6). 


het hath th (6) 


H» __ Unidirectional 
i | scratches 


Fig. 7 : Forces acting on the bubble along the surface 


A developing bubble, due to its change in momentum, 
experiences a force on its meniscus. The change in 
momentum (Fm) causes due to continual liquid is 


estimated as given in Equation 7. 
1 fa" z 
Fy = +(x‘) H, (7) 


The surface tension forces acting at top (F;,,) and bottom 
(Fp) of the bubble, are given in equation 8, respectively. 
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F,,—o0;F,, ~ ocos@ (8) 

The gravity force (F,) which acts on the boiling surface is 

a in Equation 9. 
+ 


The capillary force (F.) which acting on the bubble, as 
given in Equation 10, for capillary pressure acting inside 
the small tube of radius Re. It is assumed that 
unidirectional scratch is a capillary tube and many such 
capillary tubes may lies below the growing bubble. 


E, = ane (= 22) o cos @ (10) 


Sm 
Combining all the above forces, q^% is expressed by 
Equation 11. 
= OS 0.25 |2 
I” = hg, [og(p, — p,)1° |= + 


u + cos 8) + —— = 2a)” 


1+cos® 


(1) 
Steady-state CHF can be obtained by the Equation 12 as 
suggested by Kandlikar [11] 


1+cos 8 
TOS Ja (12) 


rr pens 
CHF. — ( 


The relation developed between steady-state CHF and Ra 
is given in Equation 13. 


Tenra = 1280.8 X R°" (13) 


The relation between q’cmrts and q°r is obtained and the 
final form of CHF model for exponential heat supply 
ranging from y=1 to y=6 for wide range of surface 
roughness R4 is given in Equation 14. 


CHF,ts = (= =] (==!) {hrap [og(e;— 


p»)]°*5} [: +=(1+cos@) + td z)” 


1+cos® 
( 14) 


Fig. 8 illustrate the comparison between experimental data 
and the predicted transient CHF values obtained from 
Equation 14 and mean absolute error (MAE) is estimated 
which found to be 4.05%. The term of y present in the 
denominator suggest that difference between steady and 
transient CHF will increase with increase in the y, 
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Fig.8 Comparison of q curs obtained by present model 
with experimental g curs 


V. CONCLUSION 


An extensive study of pool boiling on copper surface with 
wide range of surface roughness Ra under different rate of 
exponential heat supply was carried to understand the 
physical mechanism and parametric dependence of 
transient CHF. It is observed that transient CHF increases 
with increase in Ra whereas it decreases with increase in y. 
A transient CHF model is developed by including time 
constant y to account for the influence of exponential heat 
supply on the CHF. The present model predicts the 
experimental values of transient CHF with MAE=4.05%. 
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